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Abstract  _ _ _ _ 

A  model  for  introducing  the  effects  of  adiabatic  shear  banding  into  a  penetration  calculation 
was  installed  into  the  EPIC  wavecode.  These  effects  are  deemed  to  be  reduction  in  the  ratio  of 
flow  stress  to  the  value  predicted  by  the  Johnson-Cook  strength  model  and  increase  in  spall 
pressure.  A  strain-rate-  and  temperature-dependent  instability  strain  is  determined  from  small- 
amplitude  perturbation  of  constant-strain-rate  simple  shear.  Imposed  alterations  in  flow  stress 
ratio  and  spall  pressure  commence  at  the  “localization  strain,”  separated  from  the  instability  strain 
by  a  fixed  strain  increment.  The  alterations  proceed  linearly  with  increasing  effective  plastic  strain 
and  terminate  after  an  additional  fixed  strain  increment,  at  the  “failure  strain.”  The  values 
imposed  on  the  flow  stress  ratio  and  the  spall  pressure  at  the  failure  strain  are  functions  of  local 
presstne  at  the  time  step  when  localization  strain  was  reached.  The  nonzero  value  imposed  on 
the  flow  stress  ratio  in  the  case  of  positive  localization  pressure  reflects  the  phenomenon  of 
fracture  suppression  within  a  fiilly  formed  shear  band.  The  two  fixed  strain  increments  are 
evaluated  from  a  torsional  Kolsky  bar  test.  The  pre- shear-banded  spall  pressure  is  evaluated  from 
plate-on-plate  impact  data.  The  flow  stress  ratio  and  spall  pressure  at  and  beyond  the  failure 
strain  introduce  two  currently  “free”  parameters.  The  model  was  applied  to  a  set  of  problems 
involving  steel  plate  perforation  by  a  tungsten  rod,  and  reasonable  agreement  with  experiment 
was  obtainable  in  terms  of  the  final  target  hole  size  and  the  length  and  speed  of  the  tungsten 
residual. 
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1  Introduction 

Rolled  homogeneous  armor  (RHA)  is  a  medium-carbon,  martensitic,  high-strength  steel. 
U.S.  Department  of  Defense  (1991)  specifies  its  allowable  range  of  chemical  composition, 
places  broad  restrictions  on  the  heat  treatment  regimen,  and  specifies  performance  require¬ 
ments  in  terms  of  such  metrics  as  Charpy  impact  energy  and  ballistic  limit.  The  maximum 
weight  percent  of  carbon  can  be  0.30  for  plates  up  to  50.8-mm  thickness,  0.33  for  plates  be¬ 
tween  50.8-  and  101.6-mm  thickness,  and  0.35  for  plates  of  thickness  greater  than  101.6  mm. 
For  any  thickness,  the  maximum  range  in  weight-percent  carbon  is  0.10.  Some  of  the  other 
elements  present  (e.g.,  Mn,  Ni,  Cr,  Mo)  are  used  to  promote  formation  of  martensite  (Hon- 
eycombe  1981). 

In  recent  U.S.  Army  Research  Laboratory  (ARL) experiments  involving  perforation  of 
RHA  plates  by  copper  shaped  charge  jets,  cracks  in  the  remaining  RHA  target  were  lined  by 
shear  bands,  and  recovered  RHA  fragments  were  bounded  by  shear  bands  (Raftenberg  and 
Krause  1999).  Shear  banding  is  an  important  fragmentation  mechanism  in  RHA  penetration, 
yet  is  generally  not  represented  in  the  modeling.  The  result  is  a  tendency  to  underpredic  t 
target  hole  size  (Raftenberg  1996a). 

In  an  axisymmetric  penetration  calculation  with  a  Lagrangian  finite-element  wavecode, 
a  1-mm  edge  length  is  a  rough  lower  limit  to  practical  element  size.  Within  a  three-node 
triangular  element,  the  velocity  field  is  given  by  bilinear  shape  functions, 

Vr{r,z)  =  Or  +  r  -b  Cr  ^  1 
Vz\r,z)  =  +  b^r  +  C:,z  j  ' 

where  ar,br,  "  Cz  are  constants;  Vr  and  Vz  are  the  radial  and  axial  components,  respectively, 
of  velocity;  and  t  and  z  are  the  radial  and  axial  coordinates,  respectively.  The  shear  bands 
observed  in  Raftenberg  and  Krause  (1999)  had  a  thickness  of  roughly  6//m,  and  these  shape 
functions  suppress  such  fine-scaled  shear  localization.  The  goal  of  the  shear  banding  model 
is  to  introduce  into  a  1-mm  element  effects  of  shear  localization,  which  are  deemed  to  be 
reductions  in  flow  stress  and  spall  strength. 

If  shear  bands  could  be  spatially  resolved  with  a  sufficiently  fine  mesh,  perhaps  a  shear 
banding  model  would  not  be  needed  as  an  additional  ingredient.  A  flow  law  with  thermal 
softening  would  cause  local  flow  stress  reduction,  and  a  void  growth  model  with  thermal  de¬ 
pendence  would  cause  local  spall  strength  reduction.  The  former  effect,  flow  stress  reduction 
within  a  spatially  resolved  shear  band,  has  been  demonstrated  many  times  in  wavecode  sim¬ 
ulations  of  micromechanical  problems,  in  which  use  of  submicron  elements  was  practicable 
(e.g.,  Zhu  and  Batra  1991).  However,  the  practical  use  of  submicron  elements  in  macro¬ 
scopic  penetration  calculations  does  not  appear  imminent,  although  the  adaptive  meshing 
approach  (e.g.,  Camacho  and  Ortiz  1997)  is  progressing  toward  that  goal.  Moreover,  the 
spatial  resolution  of  shear  bands  would  introduce  formidable  difficulties:  extrapolations  of 
the  flow  law  and  void  model  in  terms  of  temperature,  strain,  and  strain  rate;  possible  phase 
change  modeling;  and  the  influence  of  practically  unknowable  microstructural  features  (e.g., 
individual  inclusions  and  voids)  that  can  trigger  shear  band  nucleation. 
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2  The  Model 

2.1  Overview 

The  shear  banding  model  is  shown  in  Figure  1.  This  model  is  applied  at  the  level  of  an 
individual  finite  element.  Pfaii  and  the  ratio  Y/Yjc  are  each  a  function  of  Here,  7  is  the 
imposed  flow  stress,  Pfaii  is  the  imposed  spall  pressure,  Yjc  is  the  flow  stress  as  given  by  the 
Johnson-Cook  strength  model  (Johnson  and  Cook  1983),  and  is  the  von  Mises  equivalent 
(effective)  plastic  strain,  defined  in  terms  of  e^j,  Cartesian  components  of  the  plastic  strain 

tensor,  by  _ _ 

*»  =  .  (2) 

Similarly,  Y  is  applied  to  the  von  Mises  equivalent  stress,  defined  by 


where  c'ij  are  Cartesian  components  of  the  deviatoric  Cauchy  stress  tensor  (Malvern  1969b). 
The  instability  strain,  Sinsti  is  computed  as  a  function  of  equivalent  plastic  strain  rate,  i  . 
and  temperature,  9.  User-provided  material  constants  include  the  pre-shear-banded  spall 
pressure,  Pf%,  the  s^-increment  between  initial  instability  and  the  beginning  of  localization, 
Aefoc)  fP-increment  between  the  beginning  of  localization  and  the  attainment  of  a 

fully  formed  shear  band,  The  other  user-provided  material  parameters,  b  and  Pfaii-. 

involve  the  values  attained  by  Y/Yjc  and  Pfau,  respectively,  that  characterize  the  fully 
formed  shear  band,  h  and  P^f^aii  both  functions  of  Pioc  fhe  pressure  at  the  beginning^of 
localization.  The  two  other  quantities  appearing  in  Figure  1  are  the  localization  strain, 
and  the  failure  strain,  ^fau-,  defined  by 

4oc  =  +  W 

The  features  of  the  model  are  described  and  motivated  in  the  remainder  of  §2. 


2.2  Instability  Strain 

Bai  (1982)  studied  analytically  the  growth  of  an  infinitesimal  perturbation  to  a  constant- 
strain-rate  simple  shearing  motion  of  a  viscoplastic  solid.  He  found  an  approximate  condition 
for  perturbation  growth  to  be  strain  in  excess  of  the  level  corresponding  to  maximum  stress 
on  the  applicable  constant-strain-rate  adiabatic  stress-strain  curve  (Figure  2).  This  critical 
strain  is  identified  with  shear  banding  model. 

Evaluation  of  s^^st  for  a  specific  flow  law  proceeds  as  follows.  Consider  a  flow  law  of  the 

form 

Y  =  Y{£^,£p,e)  .  (6) 

Differential  dY  along  a  general  thermodynamic  path  satisfies 


(7) 


2 


If,  during  time  increment  dt,  the  process  is  adiabatic  and  involves  nonzero  plastic  strain 
increments  (loading),  and  elastic  strain  increments  are  negligible,  then  d9  and  de^  are  related 
by 

pcde  =  0Yde^.  (8) 

P  is  the  fraction  of  plastic  work  that  is  converted  to  temperature  rise,  equal  to  about  0.9  for 
steel  (Taylor  and  Quinney  1934).^  For  such  an  adiabatic  process. 


(dY) 


'dY  pY  dY' 


dY 


adiabatic  — 


+  —  —  ]de^  +  -^de^. 


(9) 


adiabatic  pQ  QQ 


(10) 


^deP  '  pc  dO 

If  the  adiabatic  process  also  involves  a  constant  strain  rate  (de^  =  0),  then 

f-) 

The  maximum  stress  on  this  constant-strain-rate  adiabatic  stress-strain  curve  occurs  when 

de^  pc  dO 

In  the  case  of  the  Johnson-Cook  flow  law  (Johnson  and  Cook  1983), 


(12) 


where  A,  B,  C,  M,  N,  and  9m  are  material  constants,  9m  has  the  physical  interpretation  of 
melt  temperature,  and  9r  is  the  room  temperature.  Evaluation  of  dYjc/de^  and  dYjc/d9 
and  substitution  into  equation  (11)  yields  the  nonlinear  algebraic  equation  for 


if  ^  M 

- 1 

1 

[a  +  B{ePf 

1  ^  (l.O  s-i)] 

- 1 

1 

S 

1 

t-H 

_ 1 

\Bj 


0, 


(13) 


where 


F{p,i^,9)  = 


pcN{9m  -  9r) 


/3BM[l+Cln(j^)]  (j^) 


M-l 


(14) 


Equation  (13)  is  solved  for  e^^st  by  numerical  iteration  at  every  time  step  in  each  finite 
element  that  has  not  already  satisfied  the  condition.  Note  that  determination  of  €\nst  makes 
use  of  material  parameters  already  introduced  by  the  flow  law  and  requires  no  additional 
parameters.  Corresponding  to  the  4340  steel  values  for  A,  B,  C,  N,  M,  p,  c,  and  9m  published 
in  Johnson  and  Cook  (1983),  values  for  efnst  tend  to  range  between  0.4  and  0.6,  depending 
on  strain  rate  and  temperature. 

This  technique  for  evaluating  raises  some  issues;  The  conditions  of  constant  strain 
rate  and  simple  shear  will  not  generally  apply  in  a  finite  element  throughout  a  penetration 
calculation.  A  finite  amplitude  perturbation  could  cause  instability  at  a  smaller  strain  (see 
Batra  1987).  Pressure  effects  are  neglected  in  the  Johnson-Cook  strength  model.  The 
application  of  von  Mises  plasticity  in  the  definitions  of  and  d  introduces  the  assumptions 
of  initial  material  isotropy  and  subsequent  isotropic  hardening  (Malvern  1969a). 


^Note  that  temperature  contributions  from  elastic  strain  are  not  assumed  negligible  throughout  the  prob¬ 
lem,  only  during  the  increment  of  instability.  Thus,  0  remains  the  true  temperature,  with  contributions  from 
both  plastic  work  and  elastic  strain  energy. 
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2.3  Localization  Strain  and  Failure  Strain 

Marchand  and  Duffy  (1988)  performed  torsional  Kolsky  bar  tests  on  thin- walled  cylinders 
of  HY-100  steel.  Their  dynamic  torsional  stress-strain  curve  from  a  test  performed  at  room 
temperature  and  at  a  shear-strain  rate  of  1600  s"^  is  shown  in  Figure  3.  On  this  curve, 
maximum  shear  stress  corresponds  to  a  shear  strain  of  about  0.27,  or  an  effective  strain  of 
0.31,  somewhat  smaller  than  the  range  of  values  typically  obtained  for  e^nst  the  solution 
of  equation  (13)  with  the  4340  material  constants  published  in  Johnson  and  Cook  (1983). 
This  discrepancy  may  be  attributable  to  material  differences  between  HY-100  and  4340, 
and  also  to  finite  amplitude  perturbations  in  the  Marchand  and  Duffy  experiments.  Such 
finite  amplitude  perturbations  may  have  originated  from  material  microstructure,  but  also 
spuriously  from  geometric  imperfections  in  the  walls  of  the  torsion  specimen.  The  possible 
influence  of  geometric  imperfections  was  noted  in  Marchand  and  Duffy  (1988),  but  no  mea¬ 
surement  of  surface  variations  was  provided.  If  such  measurements  of  surface  imperfections 
were  available,  the  analysis  by  Molinari  and  Clifton  (1987)  could  be  used  to  estimate  the 
degree  to  which  shifted  to  smaller  values. 

In  Figure  3,  the  localization  shear  strain,  at  which  shear  stress  begins  to  decrease  signifi¬ 
cantly,  is  roughly  0.38,  which  corresponds  to  an  equivalent  strain  of  0.44.  Hence,  an  estimate 
for  Ael,  is  0.44  -  0.31  =  0.13. 

In  Figure  3,  once  localization  is  exceeded,  the  axis  label  of  “nominal  shear  strain”  is 
pertinent.  A  grid  on  the  outer  wall  of  the  specimen  was  used  to  measure  strain  near  the 
shear  band,  but  a  more  local  strain  measurement  within  the  band  was  not  obtained.  The 
figure  shows  that  by  a  nominal  shear  strain  of  about  0.57,  the  stress  is  near  zero  and  has 
effectively  ceased  to  decrease.  The  corresponding  equivalent  strain  is  0.66,  which  leads  to 
an  estimate  for  Ae^j^n  of  0.66  -  0.44  =  0.22.  Since  the  strain  measurement  was  not  ideally 
local,  this  Ae^f^n  estimate  is  a  lower  bound. 

Further  torsional  Kolsky  bar  testing  may  reveal  Aef^^  to  t)e  significantly  de¬ 

pendent  on  £,  0,  and  P.  In  the  absence  of  such  data,  they  are  here  treated  as  material 
constants.  However,  since  is  a  function  of  e  and  6,  and  with  i 

and  9  according  to  equations  (4)  and  (5). 

In  a  followup  to  Marchand  and  Duffy  (1988),  Cho,  Chi,  and  Duffy  (1990)  studied  the 
mechanism  of  strength  reduction  to  zero  within  a  shear  band  in  three  different  steels,  includ¬ 
ing  AISI  4340  with  RHC  44,  which  is  reasonably  similar  to  RHA.  (RHA  has  less  carbon  and 
an  RHC  of  about  30.)  They  performed  fractography  on  cracks  within  shear  bands  generated 
in  torsional  Kolsky  bar  specimens.  In  the  case  of  RC-44  4340,  cracking  within  the  shear 
bands  was  associated  with  the  coalescence  of  ductile  voids  that  nucleated  at  debonding  sites 
between  the  steel  matrix  and  carbide  inclusions.  The  debonding  and  subsequent  growth  of 
the  voids  were  driven  by  the  massive  shear  flow  within  the  band.  This  indicates  that  in  a 
torsion  test,  in  which  a  state  of  simple  shear  is  closely  approximated,  strength  loss  within 
a  shear  band  has  two  contributions:  (1)  thermal  softening  and  (2)  microcracking  associated 
with  ductile  voids. 
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2.4  Pressure-Dependent  Residual  Strength 

In  a  penetration  situation,  a  compressive  shock  wave  is  generated  in  both  target  and  pen- 
etrator  upon  impact.  If,  for  example,  a  tungsten-heavy-alloy  rod  impacts  an  RHA  target 
at  1.5  km/s,  the  shock  pressure  in  the  target  is  roughly  64  GPa,^  or  an  order  of  magnitude 
larger  than  the  RHA  flow  stress  of  1  to  2  GPa  and  spall  stress  of  about  6  GPa  (Bless  1981). 
This  large  superimposed  compression  is  a  feature  not  present  in  the  torsion  tests  of  Marchand 
and  Duffy  (1988)  and  Cho,  Chi,  and  Duffy  (1990). 

The  modified  torsional  Kolsky  bar  tests  of  Chichili  and  Ramesh  (Chichili  1997;  Chichili, 
Ramesh,  and  Hemker  1998;  Chichili  and  Ramesh  1999)  shed  light  on  effects  of  pressure  upon 
shear  banding.  In  these  tests,  which  involved  alpha-titanium,  an  axial  compression  wave  was 
applied  prior  to  arrival  of  the  torsional  wave.  The  specimens  were  relatively  thick- walled,  as 
needed  to  prevent  buckling  under  the  compressive  loading,  and  contained  an  axisymmetric 
notch  on  their  outer  wall.  Near  the  notch  tip,  a  finite  element  analysis  showed  the  stress  state 
to  be  approximated  by  simple  shear  plus  hydrostatic  stress.  Chichili  and  Ramesh  found  that 
without  superimposed  pressure,  the  torsion  produced  intra-shear-band  cracking  associated 
with  void  nucleation  and  growth  along  grain  boundaries.  The  application  of  compressive 
hydrostatic  stresses  on  the  order  of  the  flow  stress  were  able  to  suppress  this  cracking  within 
the  shear  band,  thereby  allowing  for  some  residual  strength.  This  ability  to  suppress  total 
strength  loss  by  means  of  pressure  reinforces  the  conclusion  arrived  at  based  on  Cho,  Chi, 
and  Duffy  (1990),  namely,  that  the  flow  stress  reduction  observed  by  Marchand  and  Duffy 
and  displayed  in  Figure  3  had  two  distinct  contributions:  thermal  softening  and  ductile 
voids. 

The  feature  of  pressure-dependent  residual  strength  is  represented  in  the  model  by  the 
function  5(P/oc),  given  by 


KPloc) 


0  ;Poc<0 

bo  e  [0, 1]  ;  Pioc  >  0 


(15) 


The  user-provided  material  constant  bo  introduces  the  assumption  that  for  strain 

hardening  in  the  fully  formed  shear  band  follows  the  Johnson-Cook  strength  model,  but  with 
amplitudes  of  deviatoric  stress  components  reduced  by  the  factor  bo- 


2.5  Spall  Strength  Reduction 

When  a  projectile  impacts  a  target,  stress  waves  are  generated  in  both  materials.  The  target 
stress  wave  reflects  upon  impact  with  the  rear  surface  to  form  a  tensile  wave.  Simultaneously, 
unloading  waves  form  along  the  penetration-hole  boundary  as  new  free  surface  is  generated. 
These  unloading  waves  interact  with  the  reflected  wave  to  form  regions  of  large  tensile 
stress,  which  can  exceed  the  material’s  spall  strength.  (Recall  that,  in  the  case  of  a  tungsten 
rod  striking  RHA,  the  initial  shock  pressure  exceeds  the  RHA  spall  stress  by  an  order  of 
magnitude  in  terms  of  absolute  value.)  Hence,  the  phenomenon  of  spallation,  or  dynamic 
failure  attributeable  to  tensile  waves,  can  occur  in  the  target. 

^This  estimate  was  obtained  from  the  relationship  Pshock  ~  Po^shock'^st  where  Pshock  is  shock  pressure, 
po  is  undeformed  RHA  density,  or  about  7800  kg/m^  (Johnson  and  Cook  1983),  Vs^ock  is  longitudinal  wave 
speed  in  RHA,  or  about  5.5  km/s,  and  u*  is  the  striking  speed. 
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Bless  (1981)  measured  spall  stress,  crfau,  in  RHA  by  plate-on-plate  impact  tests  and 
reported  the  value  of  6.0  GPa.  (Note  that  afau  is  an  axial  stress,  not  a  von  Mises  equivalent 
stress.)  Since  the  growth  of  a  spherical  ductile  void  is  generally  modeled  to  be  driven  by 
pressure  (e.g.,  Rajendran,  Dietenberger,  and  Grove  1989),  in  the  shear  banding  model  the 
spallation  criterion  is  in  terms  of  pressure.  For  the  uniaxial  strain  condition  of  plate-on-plate 
impact  tests,  pressure  and  axial  stress  are  related  by  (see  Raftenberg  1996b) 

4Y 

<yjail  +  Pfail  = 

Table  1  presents  spall  pressure,  P/aUt  for  several  reasonable  values  of  flow  stress.  In  the 
model,  these  values  are  applied  to  Pj%,  the  spall  pressure  of  pre-shear-banded  material. 

Data  on  the  effects  of  pulse  duration  on  spall  stress  are  not  available  for  RHA,  so  PjJ;  is 
taken  to  be  a  negative-valued  material  constant,  independent  of  loading  rate. 

Table  1:  Estimates  of  for  RHA  Corresponding  to  Reasonable  Values  of  Y 


y 

P fail 

(GPa) 

(GPa) 

1.0 

-4.7 

1.5 

-4.0 

2.0 

-3.3 

2.5 

-2.7 

There  is  compelling  evidence  that  shear  banding  locally  increases  (or,  in  terms  of  absolute 
value,  reduces)  spall  pressure.  Irwin  (1972)  observed  in  penetrators  composed  of  a  U-2Mo 
alloy  that  ductile  voids  formed  selectively  within  shear  bands  (Figure  4).  The  equiaxed 
nature  of  the  voids  in  this  figure  strongly  suggests  that  they  grew  under  the  influence  of 
hydrostatic  tension,  or  negative  pressure.  The  tensile  wave  was  evidently  able  to  grow  voids 
within  the  shear  band,  but  not  in  the  adjacent  material  outside  the  band.  This  phenomenon 
of  spall  pressure  increase  within  a  shear  band  is  introduced  into  the  model  by  the  function 
Pf^^iiPioc),  which  satisfies  the  relation 

vp,„.  (17) 


There  are  presumably  two  mechanisms  contributing  to  spall  pressure  increase  within 
shear  bands,  in  close  analogy  with  the  two  proposed  contributions  to  flow  stress  reduction. 
First,  spall  pressure  is  increased  by  any  void  nucleation  and  growth  that  occurred  during 
shear  band  formation.  Second,  thermal  softening  lowers  resistance  to  void  growth,  thereby 
increasing  spall  pressure.  The  first  contribution  would  again  presumably  be  eliminated  by  the 
presence  of  a  sufficiently  large  positive  pressure  during  the  process  of  shear  band  formation. 
Hence,  a  distinction  is  introduced  between  and  Pj^z"  ,  two  user-provided  material 

constants  defined  by 

.  p,  <  n 

fail  /,gN 

5(s6+)  "  "  ’ 


fail 


;  Pioc  >  0 
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where 


Pfl  <  -P/2r’  <  -p):?’  <  0  •  (19) 

That  is,  if  shear  banding  formed  under  compression,  spall  pressure  is  increased  (or  decreased 
in  terms  of  absolute  value)  from  Pf%  by  a  lesser  or  equal  amount  than  if  shear  banding 
formed  under  tension. 

Once  the  condition  P  <  Pfaii{e^ ,  Pioc)  is  satisfied,  that  element  can  no  longer  support 
deviatoric  stresses  or  hydrostatic  tension. 


3  An  Application 

3.1  Description  of  Problems 

In  1993-94,  ARL  performed  a  series  of  tests  in  each  of  which  an  RHA  plate  of  either  50.8-  or 
76.2-mm  thickness  was  perforated  by  a  small-L^/Ds  right  circular  cylinder  composed  of  91W' 
6Ni-3Co  tungsten  heavy  alloy  (WHA),  which  impacted  at  normal  incidence  (Raftenberg  and 
Kennedy  1995).  The  six  problems  that  were  studied  are  described  in  Table  2;  Vg,  Lg,  and 
Dg  are  the  striking  speed,  length,  and  diameter,  respectively,  of  the  WHA  rod,  and  d  is  the 
thickness  of  the  RHA  target,  a  circular  plate  with  a  1-m  diameter.  The  rear  of  the  rod  was 
surrounded  by  a  stabilizing  7075-T651  aluminum  drag  flare  in  the  form  of  a  truncated  cone. 
A  total  of  10  experiments  were  performed  (Table  3).  In  each,  the  impact  speed  was  within 
100  m/s  of  the  desired  value,  and  the  impact  obliquity,  7^,  was  3.00  degrees  or  less.  The 
experiments  are  more  thoroughly  described  in  Raftenberg  and  Kennedy  (1995). 


Table  2:  Six  Problems  Deflned 


Problem 

Vs 

(mm//zs) 

Ls 

(mm) 

Ds 

(mm) 

d 

(mm) 

1 

1.52 

112.0 

20.9 

50.8 

2 

1.90 

92.0 

18.3 

50.8 

3 

2.30 

65.0 

16.3 

50.8 

4 

1.52 

112.0 

20.9 

76.2 

5 

1.90 

92.0 

18.3 

76.2 

6 

2.30 

65.0 

16.3 

76.2 

3.2  Description  of  the  Computations 

The  six  problems  were  simulated  with  the  EPIC  lagrangian  wavecode  (Johnson  and  Stryk 
1992),  into  which  the  shear  banding  model  had  been  installed.  All  finite  elements  were 
axisymmetric  three-node  triangles  arranged  in  rectangular  groups  of  four  crossed-triangles. 
There  were  5  such  rectangles  across  the  radius  of  the  rod  and  50  along  the  length.  For  the 
50.8-mm-thick  targets,  there  were  25  such  rectangles  across  the  thickness  and  200  across  the 
radius.  For  the  76.2-mm-thick  targets,  there  were  38  such  rectangles  across  the  thickness 
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Table  3:  Ten  Experiments  Identified 


Test 

Vs 

Ls 

Ds 

Is 

d 

Problem 

No. 

(mm///s) 

(mm) 

(mm) 

(deg) 

(mm) 

1086 

1.57 

112.0 

20.9 

0.40 

50.8 

1 

1071 

1.84 

92.0 

18.3 

3.00 

50.8 

2 

1089 

1.93 

92.0 

18.3 

2.54 

50.8 

2 

1069 

2.30 

65.0 

16.3 

2.33' 

50.8 

3 

1087 

1.59 

112.0 

20.9 

2.24 

76.2 

4 

1076 

1.62 

112.0 

20.9 

0.90 

76.2 

4 

1090 

1.91 

92.0 

18.3 

1.18 

76.2 

5 

1075 

1.99 

92.0 

18.3 

1.56 

76.2 

5 

1092 

2.26 

65.0 

16.3 

2.22 

76.2 

6 

and  200  across  the  radius.  Eroding  slidelines  (Johnson  and  Stryk  1996)  were  located  at  the 
interfaces  between  rod  and  target,  rod  and  flare,  and  flare  and  target;  an  erosion  strain  of 
1.5  and  a  Coulombic  friction  coefficient  of  zero  were  used  throughout. 

The  shear  banding  model  was  applied  to  the  RHA  target  material  in  all  calculations,  with 
the  parameters  Pj^u  set  to  —3.0  GPa  (based  on  Table  1),  to  0.13,  and  to  0.22 

(the  last  two  based  on  Figure  3).  and  PjJr’  were  always  equated  (hence,  the  symbol 

is  used  as  replacement),  and  their  shared  value  was' varied  in  the  range  between  -3.0 
—0.2  GPa.  bo  was  varied  between  0.0,  0.5,  and  1.0.  Note  that  bo  —  0.0  is  the  special 
case  when  the  fully  formed  shear  band  has  no  residual  strength,  bo  =  1.0  is  the  case  when 
Yjc  is  applied  to  the  fully  formed  shear  band  (except  in  the  rare  occasion  that  Poe  <  0),  so 
that  the  shear  banding  model  degenerates  to  spall  pressure  increase  only. 

For  the  WHA  rod  material  and  the  aluminum  of  the  flare,  the  Johnson-Cook  fracture 
model  (Johnson  and  Cook  1985)  was  used  to  represent  damage.  This  model  introduces  the 
seven  material  constants  Di,  D2,  Dz,  D^,  D5,  a  spall-  These  were  assigned  the 

values  in  Table  4,  which  were  obtained  from  Johnson  (1997)  for  7039  A1  and  90W-7Ni-3Fe. 

The  Johnson-Cook  strength  model  was  used  to  evaluate  flow  stress,  Yjc,  for  RHA,  WHA, 
and  Al.  The  von  Mises  plasticity  algorithm  imposed  isotropic  hardening.  Elastic  deviatoric 
strains  were  related  to  deviatoric  stresses  by  the  elastic  shear  modulus,  G,  assumed  to  be 
a  material  constant.  G  and  the  Johnson-Cook-strength  constants  A,  B,  C,  M,  N,  and 
were  assigned  the  values  from  Johnson  and  Cook  (1983)  that  are  listed  in  Table  4.  The 
values  given  in  Johnson  and  Cook  (1983)  for  Rockwell  C-30,  austenitized,  quenched,  and 
tempered  4340  steel  were  assigned  to  RHA.  Those  given  for  90W-7Ni— 3Fe  were  assigned  to 
the  91W-6Ni-3Co  WHA,  and  those  for  7039  Al  were  assigned  to  the  7075-T651  Al. 

Dilatation  of  the  three  metals  was  governed  by  the  Mie-Griineison  equation  of  state,  with 
the  shock-Hugoniot  pressure,  Ph,  related  to  compression,  p,  by  the  cubic 

Pniy)  =  ,  (20) 


where 


H  = 


P  -  Po 


Po 


(21) 
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Table  4:  Material-Constant  Values 


Material 

Constant 

A1 

WHA 

RHA 

Po  (kg/m'") 

2768. 

17000. 

7830. 

Ki  (GPa) 

76.74 

246.1 

163.9 

K2  (GPa) 

128.3 

391.9 

294.4 

Kz  (GPa) 

125.1 

820.8 

500.0 

r 

2.00 

1.43 

1.16 

G  (GPa) 

26.20 

124.11 

77.50 

A  (GPa) 

0.3365 

1.5058 

0.7922 

B  (GPa) 

0.3427 

0.1765 

0.5095 

C 

0.010 

0.016 

0.014 

M 

0.41 

0.12 

0.26 

N 

1.00 

1.00 

1.03 

0riK) 

294.3 

294.3 

294.3 

0m  (K) 

877.6 

1723.2 

1793.2 

C  (J/kg-K) 

875.6 

134.5 

477.8 

Di 

0.14 

0.0 

— 

T>2 

0.14 

0.33 

— 

Dz 

-1.5 

-1.5 

— 

D, 

0.018 

0.000 

— 

Dz 

0.00 

0.00 

— 

eD 

^TTitn 

0.060 

0.022 

— 

(^spaii  (GPa) 

4.62 

6.76 

— 

(GP^) 

— 

— 

-3.00 

— 

— 

0.13 

AeLi 

— 

— 

0.22 
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p  is  the  current  density,  and  po  is  the  undeformed  density.  Ki,  K2,  and  are  material 
constants,  as  is  the  Griineison  coefficient,  T.  Values  assigned  to  these  constants  and  to  Po 
are  listed  in  Table  4.  For  aluminum,  Po  was  obtained  from  Johnson  and  Cook  (1983)  for 
2024-T351  Al;  and  Ki,  K2,  Ks,  and  F  were  obtained  from  Kohn  (1969)  for  2024  Al.  For 
RHA,  Po  was  obtained  from  Johnson  and  Cook  (1983)  for  4340  steel.  Ki,  K2,  K3,  and  F 
were  obtained  from  Kohn  (1969)  for  304  stainless  steel.  For  WHA,  po  was  obtained  from 
Johnson  and  Cook  (1983)  for  a  90W-7Ni-3Fe  alloy.  F  was  assigned  the  value  in  Kohn  (1969) 
for  pure  tungsten,  while  iFi,  iF2,  and  Kz  were  evaluated  from  the  linear  shock  speed-particle 
speed  curve  reported  in  Hauver  (1980)  for  90W-7Ni-3Fe. 


3.3  Computational  and  Experimental  Results  Compared 

After  each  experiment,  the  target  plate  was  sectioned  along  an  approximate  perforation- hole 
diameter  to  reveal  hole  morphology  and  locations  of  cracks  and  cavities.  Figure  5  shows  the 
sectioned  plate  from  Test  1071.  Note  the  cracks  that  enter  the  plate  from  the  hole  boundary. 
Also  note  the  “spall  ring  shelf”  that  surrounds  the  hole  boundary  at  the  exit  surface.  This 
feature  was  created  by  the  separation  from  the  plate  of  the  largest  fragments  that  were 
observed.  For  each  plate,  the  hole  boundary,  spall  ring  shelf,  cracks  emanating  from  the  hole 
boundary,  and  internal  cavities  were  digitized.  These  digitizations  in  the  case  of  Test  1086  are 
shown  in  Figure  6.  From  the  hole  boundary  digitization,  the  through-thickness- averaged  hole 
diameter,  D,  was  computed  by  integrating  from  the  initial  elevation  of  the  entrance  surface, 
z  =  0,  to  that  of  the  exit  surface,  z  =  -50.8  mm  or  -76.2  mm.  Note  that  this  integration  did 
not  include  contributions  from  the  spall  ring  shelf  region  since  the  lip  of  the  overhang  crossed 
the  initial  elevation  of  the  exit  surface.  An  uncertainty  in  the  experimental  D  value  of  about 
±0.3  mm  was  introduced  by  the  boundary  digitization  procedure.  Additional  uncertainties 
introduced  by  deviations  from  the  intended  impact  speed  and  small  unintended  yaw  angles 
are  difficult  to  quantify.  Flash  radiography  was  used  to  measure  the  length,  L^es,  and  speed, 
Vres,  of  the  residual  rod  shortly  after  its  emergence  from  the  target  exit  surface. 

All  EPIC  simulations  were  run  for  2.5  ms  after  initial  impact  at  time  t  =  0.  The  quantites 
D,  Lres,  and  Vres  were  compu^d  from  each  EPIC  simulation  at  2.5  ms. 

Figures  7  through  12  plot  D  vs.  P)®*]  with  bo  a  parameter  for  each  of  the  six  problems.  In 
each  of  these  figures,  the  dashed  line(s)  indicate  experimental  measurements.  The  computa¬ 
tional  point  =  -3.0  GPa,  bo  -  1-0  is  essentially  the  solution  without  the  shear  banding- 
model.  In  each  of  Problems  1  through  5  (no  D  measurement  was  obtained  for  Problem  6), 
this  point  is  significantly  smaller  than  the  measured  values.  In  Problem  1,  for  instance,  the 
computed  D  of  38.8  mm  is  11.5%  smaller  than  the  exp^imental  value  of  43.6  mm,  or  roughly 
23%  smaller  in  terms  of  hole  volume.  Figure  13  plots  D  vs.  u,  from  the  six  problems  for  this 
case  of  Pffi  =  -3.0  GPa,  bo  =  1.0. 

In  Figures  7  through  12,  for  a  given  P)^3  <  -l.OGPa  the  computed  D  results  exhibit  the 
anticipated  trend  of  increasing  as  bo  (the  residual  Y/Yjc  value  wiAin  the  fully  formed  shear 
band)  decreases  through  the  values  1.0,  0.5,  0.2,  and  0.0.  The  D  results  for  bo  =  0.0  and 
bo  =  0.2  are  relatively  insensitive  to  Pj^ii  throughout  the  range  Pjf^i  ^  0.0.  In  Problems  2 
and  4,  for  each  of  which  there  are  two  experimental  measurements  for  D,  bo  =  0.0  results  are 
close  to  the  larger  measurement,  and  bo  =  0.2  results  are  close  to  the  smaller  measurement.  In 
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Problem  5,  for  which  there  are  also  two  experimental  measurements  for  D,  bo  =  0.0  results  lie 
between  the  two  measurements  and  the  bo  =  0.2  results  are  slightly  smaller  than  the  smaller 
measurement.  The  scatter  of  2  to  3  mm  between  the  two  experimental  measurements  in 
these  figures  can  perhaps  be  attributed  to  variations  in  impact  speed  and  yaw  angle.  In 
Problem  1,  the  bo  =  0.0  and  bo  =  0.2  results  lie  above  the  single  measurement,  while  m 
Problem  3  they  lie  below  the  single  measurement.  Throughout  Problems  1  through  6,  D 
results  corresponding  to  bo  =  0.5  and  bo  =  1.0  are  rather  insensitive  to  throughout 

the  range  <  -1.0  GPa,  but  then  increase  with  increasing  for  P{2  >  -1.0  GPa. 

Throughout  Problems  1  through  5,  at  P}2  =  -0.2  GPa  the  D  results  for  bo  =  0.0  and 
=  0.5  come  reasonably  close  to  one  or  more  measured  values. 

Figure  14  plots  Lres  vs.  P)®^;  with  bo  a  parameter  for  Problem  2.  For  a  given  problem, 
one  or  two  measurements  were  obtained  for  the  experimental  residual  rod  length, 
while  24  values  were  obtained  for  computational  residual  rod  length,  corresponding 

to  permutations  of  Pf^au  ~  —3.0,  —2.0,  —1.5,  —1.0,  —0.5,  —0.2  GPa  and  bo  =  0.0, 0.2, 0.5, 1.0. 

is  defined  to  be  the  maximum  discrepancy  between  an  value  and  a  specific 

gQ 

=  max  ^  \LZ7  -  4SI  •  (22) 

(«.».) 

l\comp-exp\  values  from  the  nine  experiments  for  which  an  x-ray  flash  allowed  for  the  deter¬ 
mination  of  are  listed  in  Table  5.  is  defined  to  be  an  experimentally  measured 

net  rod  shortening, 

AL^^p  =  Ls-LfTs-  (23) 

In  Table  5  is  normalized  by  AL^^^,  and  we  see  that  the  computational  “error”  in 

residual  rod  length  varies  between  4.1  and  18.7%  of  the  net  shortening.  Hence,  reasonabh' 
accurate  results  for  rod  shortening  were  obtained  from  all  EPIC  runs,  including  those  for 
which  the  shear  banding  model  was  effectively  inoperative  (Pjail  =  ~3.0  GPa,  bo  =  1.0). 


Table  5:  Comparison  of  Computational  and  Experimental  Lres  Results 


Problem 

Ls 

AL^^p 

T  \comp—exp\ 
■^res 

T-  \comp—exp\ 
^res 

AL^^P 

(mm) 

(mm) 

(mm) 

(mm) 

(%) 

1 

112.0 

67.2 

44.8 

2.5 

5.6 

2 

92.0 

48.8 

43.2 

4.7 

10.9 

46.0 

46.0 

3.3 

7.2 

3 

65.0 

23.4 

41.6 

6.1 

14.7 

4 

112.0 

49.3 

62.7 

2.6 

4.1 

52.6 

59.4 

5.9 

9.9 

5 

92.0 

33.1 

58.9 

5.4 

9.2 

34.0 

58.0 

6.3 

10.9 

6 

65.0 

16.3 

48.7 

9.1 

18.7 

Figure  15  plots  Vres  vs.  Pj^-J  for  Problem  2.  As  shown  in  Tables  2  and  3,  for  some 
experiments  there  was  considerable  discrepancy  between  the  achieved  striking  speed, 


11 


(measured  with  a  streak  camera)  and  the  intended  value.  From  each  of  the  seven  exper¬ 
iments  for  which  two  x-ray  flashes  allowed  for  a  determination  of  residual  speed,  an 
experimental  net  rod  deceleration,  defined  by 


^  ^exp_^exp 


(24) 


was  determined  (Table  6).  For  each  problem,  the  striking  speed  in  the  EPIC  simulations, 
ycomp^  of  course  the  intended  value,  and  corresponding  to  each  of  the  24  permutations 
of  {plfl  bo),  a  residual  speed,  was  computated.  The  minimum  and  maximum  results 

for  rod  deceleration,  and  respectively,  defined  by 

=  max  {vr^  -  v^eTn  ,  (25) 

,  min  (vr^-vreD,  (26) 

are  also  listed  in  Table  6.  Comparison  of  A?;^j„^  and  Au^™^  with  An®  ^  in  Problems  3 
and  4  reveals  that  computational  results  for  rod  deceleration  were  significantly  smaller  than 
experimental  results  for  these  two  problems.  In  Problem  1,  computational  decelerations  agrc'o 
with  one  of  the  two  experiments  but  are  smaller  than  results  from  the  other.  In  Problem  5, 
computational  decelerations  are  bounded  by  results  from  the  two  experiments. 


Table  6:  Comparison  of  Computational  and  Experimental  Vres  Results 


Problem 

vf^ 

yexp 

'^res 

comp 

A^.comp 
— *  ^max 

(mm/^s) 

(mm///s) 

(mm//xs) 

(mm/fjis) 

(mm//is) 

(mm/)us) 

1 

1.57 

— 

— 

1.52 

0.12 

0.14 

2 

1.84 

1.63 

0.21 

1.90 

0.15 

0.19 

1.93 

1.74 

0.19 

3 

2.30 

1.94 

0.36 

2.30 

0.22 

0.32 

4 

1.59 

1.29 

0.30 

1.52 

0.22 

0.26 

1.62 

1.35 

0.27 

5 

1.91 

1.58 

0.33 

1.90 

0.25 

0.32 

1.99 

1.77 

0.22 

6 

2.26 

— 

— 

2.30 

0.45 

0.67 

Figures  6  and  16  through  20  show  a  time  sequence  of  mesh  plots  from  Problem  1  with 
=  -0.2  GPa  and  bo  =  0.5.  The  finite  element  meshes  are  superimposed  over  the 
digitizations  of  experimental  hole  profile  and  crack  and  cavity  locations.  The  Computational 
Damage”  legend  indicates  the  status  of  colored  finite  elements  with  regard  to  the  shear 
banding  model.  Table  7  interprets  legend  entries  in  terms  of  the  model.  At  25  after 
impact  (Figure  16),  there  are  small  regions  of  spalled  RHA  elements  lining  the  perforation 
hole.  Among  the  RHA  elements  that  have  not  yet  spalled,  those  closest  to  the  hole  satisfy 
the  condition  6%^  <  e^.  Adjacent  to  these  is  a  band  of  elements  that  satisfy  <e^  <  e)ou^ 
and  these  are  followed  by  a  band  that  satisfy  4,,  <  <  eL-  At  50  /iS  (Figure  17),  there 
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is  a  small  region  of  RHA  near  the  exit  surface  that  has  spalled,  satisfying  the  condition 
<  ^inst  ;  ^  Pfaii-  Figure  20  that  reasonable  agreement  has  been  achieved 

in  terms  of  P,  but  the  spall  ring  shelf  is  missing  from  the  calculated  result.  Also,  colored 
elements  that  have  satisfied  the  instability  condition  do  not  extend  as  far  radially  into  the 
targets  as  do  the  experimental  cracks,  each  one  of  which  presumably  runs  along  a  shear  band 
(Raftenberg  and  Krause,  1999). 


Table  7:  Code  for  Damage  Legend  on  Mesh  Plots 


Legend  Entry 
NO  INSTAB;  SPALL 
INSTAB  BUT  NO  LOG;  NO  SPALL 
INSTAB  BUT  NO  LOG;  SPALL 
LOG  BUT  NO  FAILURE;  NO  SPALL 

LOG  BUT  NO  FAILURE;  SPALL 
FAILURE;  NO  SPALL 
FAILURE;  SPALL 


Interpretation 


q  < 

^tnst 

7 

P 

< 

p{^) 

^fail 

^inst 

< 

<  ; 

P 

> 

p(o) 

'^fail 

^znst 

< 

<  ; 

P 

< 

p{^) 

^fail 

^loc  - 

P 

> 

p(o) 

+ 

rp(s6)  _ 
l-^fail 

-  p^°'^  1 

^loc  ■ 

P 

< 

p(o) 

+ 

\  p{s^) 
[■^fail 

-  1 
^faiU 

^fail 

< 

5 

P 

> 

^fail 

^  fail 

7 

P 

< 

p(^^) 
^ fail 

Figures  21  through  23  present  computational  and  experimental  results  from  the  six  prob¬ 
lems  when  =  -0.2  GPa,  K  =  0.5.  Comparison  of  Figure  21  with  Figure  13  shows  that 
the  shear  banding  model  has  allowed  for  improved  agreement  with  experimental  hole  size 
while  maintaining  reasonable  agreement  with  measured  values  for  Lres  and  Vres  (Figures  22 
and  23). 


4  Concluding  Remarks 

4.1  Summary  of  Results 

A  shear  banding  model  was  developed  to  impose  on  a  finite  element  the  effects  of  adiabatic 
shear.  These  are  deemed  to  be  reduction  in  the  ratio  of  flow  stress  to  the  value  predicted 
by  the  Johnson-Cook  strength  model  and  increase  in  spall  pressure.  The  instability  strain, 
,  is  identified  with  the  maximum  on  the  constant-strain-rate,  adiabatic  stress-strain  curve 
at  the  current  level  of  strain  rate.  After  a  strain  increment,  the  localization  strain 

is  reached,  at  which  point  the  reduction  in  the  flow  stress  ratio  and  the  increase  in  spall 
pressure  from  its  initial  value  of  both  commence.  This  reduction  and  increase  proceed 
as  linear  functions  of  equivalent  plastic  strain  until  a  second  strain  increment,  is 

achieved,  at  which  point  the  flow  stress  ratio  attains  its  residual  value,  b{Pioc),  and  the  spall 
pressure  attains  its  shear-banded  value  of  Pj^2  treated  as  a  material  constant, 

and  b{Pioc)  introduces  another  material  constant,  bo- 
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Afff  Ae^t  ;  and  Pf°2i  were  also  treated  as  material  constants.  The  first  two  were 
assigned  values  from  a  constant-strain-rate  torsional  Kolsky  bar  test  on  HY-100,  and  the 
third  from  a  plate-on-plate  impact  test  on  RHA.  and  bo  are  at  present  unevaluated  and 
were  treated  as  “free”  parameters. 

The  model  was  installed  into  EPIC  and  applied  to  a  set  of  six  problems  involving  perfora¬ 
tion  of  RHA  plates  by  WHA  rods.  The  model  was  found  able  to  rectify  an  important  short¬ 
coming  in  previous  calculations,  namely,  a  sufficiently  large  target  hole  size  was  achieved, 
albeit  by  means  of  “free”  parameters.  The  feature  of  the  spall  ring  shelf  at  the  target  exit 
surface  remains  elusive. 

4.2  Suggestions  for  Future  Work 

For  the  application  of  RH.A  plate  perforation,  evaluations  of  the  model  s  parameters  Asi^^. 
and  Ae^foii  can  be  refined  by  means  of  torsional  Kolsky  bar  tests  on  RHA.  The  possible 
dependency  of  these  parameters  upon  strain  rate  can  also  thereby  be  explored. 

The  currently  unevaluated  parameters  and  bo  have  a  clear  physical  significance  and 
should  in  principle  be  amenable  to  experimental  evaluation.  The  compression-torsion  Kolsky 
bar  (Chichili  1997;  Chichili  and  Ramesh  1999)  may  provide  a  means  for  evaluating  bo  during 
the  compression  phase  of  the  experiment,  when  cracking  within  the  shear  band  is  suppressed. 
However,  much  uncertainty  is  introduced  by  the  finite  element  calculation  that  is  used  to 
estimate  residual  stress  within  the  shear  band  (Ramesh  1999).  Using  this  same  apparatus, 
perhaps  by  not  trapping  the  compression  wave  a  tensile  wave  can  be  subsequently  delivered 
to  the  fully  formed  shear  band  with  some  controlled  delay  of  arrival  time,  and  evaluation  of 
Pif2i  can  thereby  be  investigated. 

Kerley  (1993)  simulated  by  means  of  the  Eulerian  CTH  code  (McGlaun  et  al.,  1990) 
an  experiment  conducted  by  Raftenberg  (Raftenberg  1994),  in  which  a  13-mm-thick  RHA 
plate  was  perforated  by  a  copper  shaped  charge  jet.  The  simulation  included  the  generally 
neglected  phenomenon  of  the  e  iron  phase  transformation  (Kerley  1993),  and  this  feature 
was  found  crucial  in  obtaining  an  accurate  hole  morphology,  particularly  at  the  regions  of 
the  entrance  and  exit  faces.  (However,  the  calculation  used  perfect  plasticity  and  contained 
no  damage  model,  and  so  was  simplistic  in  other  regards.)  The  thin  13-mm-thick  plate 
contained  no  true  spall  ring  shelf  in  the  simulations  or  experiment,  but  the  calculations 
are  nevertheless  suggestive  that  for  a  thicker  plate  inclusion  of  the  phase  transformation 
representation  may  be  what  is  needed  to  produce  the  shelf. 


14 


Figure  3:  Torsional  Stress-Strain  Curve  for  HY-100  Steel  at  Room  Temperature  and  a  Con¬ 
stant  Strain  Rate  of  1600  s“^  (reproduced  from  Marchand  and  Duffy  [1988],  Figure  12,  with 
permission  from  the  publisher). 
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SHEAR  STRESS  -  MPo 


Figure  4:  Ductile  Voids  Within  a  Shear  Band  in  a  U-2Mo  Penetrator  (reproduced  from  Irwin 
[1972]  with  permission  from  Defence  Research  Establishment  Valcartier). 
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Figure  6:  Initial  Mesh  for  Problem  1,  Superimposed  on  Digitized  Experimental  Perforation 
Hole  Contour  With  Cavities  and  Cracks  Indicated. 
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Figure  7:  Final  Through-Thickness-Averaged  Hole  Diameter  vs.  Spall  Pressure  of  Shear- 
Banded  RHA  From  Problem  1. 
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Figure  8:  Final  Through-Thickness-Averaged  Hole  Diameter  vs.  Spall  Pressure  of  Shear- 
Banded  RHA  Prom  Problem  2. 
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Figure  9:  Final  Through-Thickness-Averaged  Hole  Diameter  vs.  Spall  Pressure  of  Shear- 
Banded  RHA  From  Problem  3. 
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Figure  10:  Final  Through-Thickness-Averaged  Hole  Diameter  vs.  Spall  Pressure  of  Shear- 
Banded  RHA  From  Problem  4. 
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Figure  11:  Final  Through-Thickness- Averaged  Hole  Diameter  vs.  Spall  Pressure  of  Shear- 
Banded  RHA  From  Problem  5. 


25 


Problem  6 


(GPa) 


Figure  12:  Final  Through-Thickness-Averaged  Hole  Diameter  vs.  Spall  Pressure  of  Shear- 
Banded  RHA  From  Problem  6. 
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Figure  13:  Final  Through-Thickness- Averaged  Hole  Diameter  vs.  Striking  Speed  From  Prob¬ 
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Figure  14:  Residual  Rod  Length  vs.  Spall  Pressure  of  Shear-Banded  RHA  From  Problem  2. 
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Figure  15:  Residual  Rod  Speed  vs.  Spall  Pressure  of  Shear-Banded  RHA  From  Problem  2. 
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Figure  18:  Mesh  Plot  at  75  pis  After  Impact  From  Problem  1  with  =  “3.0  GPa  and 
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Figure  20;  Mesh  Plot  at  2500  fts  After  Impact  From  Problem  1  with  =  -3-0  GPa  and 
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Figure  21:  Final  Through-Thickness- Averaged  Hole  Diameter  vs.  Striking  Speed  From  Prob¬ 
lems  1  Through  6  for  =  —0.2  GPa  and  bo  =  0.5. 
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Figure  22:  Residual  Rod  Length  vs.  Striking  Speed  From  Problems  1  Through  6  for 
=  -0.2  GPa  and  K  =  0.5. 
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Figure  23:  Residual  Rod  Speed  vs.  Striking  Speed  From  Problems  1  Through 
pg  =  -0.2  GPa  and  h  =  0.5. 
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